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1. Introduction 

Interest in SU(N) gauge theories witli large N began with the pioneering studies of [|l|], who 
showed that in the limit N ^ oo and the gauge coupling, g — )• 0, with the 't Hooft coupling A = g^N 
fixed, one gets a non-trivial but simplified theory. Many qualitative features of hadron physics can 
be explained by appealing to this limit. 

The theory with an infinite number of colors has also been used to understand various facets 
of the phase diagram of strongly interacting matter. At very high temperatures, strongly interacting 
matter is known to exist in a deconfined, chirally symmetric state. The nature of the transition to 
this state is quite sensitive to the quark sector. For infinitely massive quarks, one has a first order 
deconfinement transition for N > 2. For two massless quarks, on the other hand, one has a second 
order, chiral symmetry restoring transition. If a small mass is given to the quarks, one expects a 
crossover, and a critical point at some nonzero baryon density. See Ref. iQ] for a detailed discussion 
of the nature of the transition. A rich phase structure has also been predicted for the large baryon 
density regime. In particular, a chirally symmetric, confined phase has been predicted, using the 't 
Hooft limit ig]. 

For phase diagrams based on arguments about the theory with an infinite number of colors to 
be relevant for the theory with three colors, it is important to explore the validity of the large N 
arguments for N=3. Analogies have also been drawn between the high temperature phase of QCD 
and the solvable, conformal ^ =4 supersymmetric SU(N — )• oo) theory Lattice studies of the 
theory with moderate values of N can provide one intermediate step in this connection, by giving 
an estimate of the size of the 1 /N^ corrections. 

In this report, we present results for a study of the thermodynamics of the SU(N) gauge theories 
at finite temperatures for N = 4,6, and combine them with results for N=3 to get an estimate of the 
applicability of the large N arguments for the theory with three colors. For N=3, we have used 
the simulation results of Ref. [^, supplementing them with new simulations where necessary (in 
particular, for the measurement of the latent heat), and analyzed them using the same techniques as 
used for the N=4 and 6 theories. For estimating the thermodynamic quantities, we use the standard 
methods The pressure, p, and A = £ — 3p, where e is the energy density, are calculated from 
the plaquette data, 

where 5P(j8, T) = P{^,T) —P{^,T = 0) is the difference in the plaquette observables between the 
finte temperature lattice and the corresponding zero temperature lattice, calculated at the coupling 
j8. To is some reference temperature. dj5 /dloga is related to the beta function. We use a nonpertur- 
bative estimate for the beta function, using the scaling of from [^. Details of the beta function 
used by us can be found in [Q]. 

In the next section we discuss the latent heat of the deconfinement transition for the theory with 
N=3,4,6 number of colors. In Sec. ^ we study the equation of state of the SU(N) gluon plasma. 
A comparison of the energy density and pressure for SU(3,4,6) gives us an idea of the size of the 
1 /N^ corrections. In the final section we discuss the physical implications of our results. 
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SU(N) gauge theories have been investigated previously on the lattice. In particular, the latent 
heat of the transition has been studied in Ref. ^, while the equation of state has been studied 
in Ref. |l^]. The focus of our work is in getting results for the continuum limit. We use a 
nonperturbative beta function, and multiple lattice spacings at each temperature, for this purpose. 
We have also employed several spatial volumes at each lattice spacing, to reach the thermodynamic 
limit. Details of our work, including the cutoff and volume dependence, can be obtained in Ref. 
[^]. The results reported here are our estimates of the infinite volume, continuum results. An 
earlier version of our study, which used the two-loop beta function, and had measurements at fewer 
temperatures, was presented in the previous year's conference [^. 



2. Latent Heat 



SU(N) gauge theories have a first order deconfinement transition for N > 3. A first order 
transition is characterized by a latent heat associated with the transition. For the infinite volume 
system, the latent heat can be defined as 

where e{Tc it 5T) are the energy densities of the confined and deconfined phases, respectively, and 
the second equality follows from the fact that pressure is continuous across the transition. 

For a finite volume system, as is necessarily used in a numerical lattice computation, one 
cannot get a separation of the phases, and so, a straightforward 5r — )• limit will not work. We use 
the following method to extract the latent heat [Q]. Since the confined and the deconfined phases 
are resolved by the Polyakov loop |L|, we identify the configurations at the transition point with 
\L\ < Lc as being in the confined phase and those with \L\ > L/, as being in the deconfined phase, 
where Lc and L/; are suitably chosen values. In order for the procedure to be meaningful and not 
too sensitive on the choice of L/, and Lc, it is important that \L\ shows a two-peak structure. Figure 
[I] (a) shows the distribution of \L\ in the transition regime of SU(4) gauge theory. Lc and L/, are 
chosen at the valley (see Ref. [§] for details). The distribution of dP{T) for the two phases defined 
by the \L\ cut is quite stable, as shown in Fig. |l](b): for values of the coupling j3 in the transition 
regime, the histograms for 5P{T) for the two phases defined by the \L\ cut essentially coincide. 

From the histograms of 5P{T) for the two phases at Tc, the latent heat can be calculated using 
Eq.Q and Eq.(|n|): 



where, as mentioned in the introduction, dp/d\oga is evaluated from the nonperturbative scale- 
setting in We find that the finite volume effects in the latent heat are rather large for SU(3), 
where one needs an aspect ratio = LT > 6 to reach the thermodynamic limit. On the other hand, 
for SU(4) and SU(6) gauge theories an aspect ratio > 3 suffices to reach the thermodynamic 
limit. This is probably related to the large correlation length of the Polyakov loop for SU(3) gauge 
theory near the transition. Our estimates for the thermodynamic limits of the latent heat for the 
theories with N = 3, 4 and 6 are shown in Table To estimate the sensitivity of the results on the 
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Figure 1: (a) Distribution of the Polyakov loop, ILI, in the transition region for SU(4) gauge theory, (b) 
Distribution of the plaquette, P{T) — P{T = 0), for the cold and hot phases, identified by configurations with 
|L| < Lc and \L\ > L^, respectively. 



Table 1: Latent heat for SU(N) gauge theories with N=3, 4 and 6. The numbers in brackets are the errors on 
the least significant digits; the first one is statistical and the second, systematic. 



N 



1.67(4)(4) 4.32(6)(6) 11.93(34)(5) 



choices of L^/,, we vary L/, — Lc by ±20%; the corresponding change in the latent heat is shown as 
a systematic error in the table. 

Our results for the latent heat are in good agreement with Ref. [j8|], within the larger uncertain- 
ties of that study. The values obtained by Ref. [Q] are somewhat higher. On going from N=3 to 4, 
the latent heat is seen to scale faster than the Casimir = — 1, the dimensionality of the adjoint 
representation. On the other hand, the scaling between N=4 and 6 is consistent with a Casimir 
scaling. Taking all the three values in the Table, one can get a good fit to the relation 

- 0.388(3) - (2-3) 



where the errors are statistical only. Eq.( |2.3| ) indicates a large correction to the leading term for the 
theory with three colors. 

3. Equation of State 

We first discuss A, which is the trace of the energy-momentum tensor and therefore is a mea- 
sure of conformal symmetry breaking. Fig. |^ shows A for SU(N) with different number of colors, 
measured at different temperatures, and normalised by c^r^. 

As the figure shows, A scales nicely with dA except very close to Tc. The size of the 1 /N^ 
correction is smaller than the statistical accuracy of our data, for T > 1.25 T^. Closer to Tc we see 
a deviation from scaling with dA- The peak of A is seen to become higher and move closer to Tc 
with increasing N. At higher temperatures, A scales approximately like T^ 1 12], for the temperature 
range investigated by us [0]. 
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Figure 2: A/c/^T"* for SU(3-6) gauge theories, plotted against Tc/T. For SU(3), we have used the plaquette 
dataofRef. 

The pressure was calculated using the so-called integral method [^], Eq.(|l.l[), which deter- 
mines p{T) in terms of pressure at some reference temperature p{To). We find that p{T) ~ 
within our errors till temperatures ~ 0.97^ . This is probably related to the fact that glueballs are 
much heavier than Tc, and therefore are not excited substantially except very close to T^. We there- 
fore evaluate p/T'^ by taking po{T < O.Sr^.) as the lower limit of the integral, with the knowledge 
that the omitted additive constant is smaller than the statistical error. 
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Figure 3: (a) The pressure and energy density for SU(N) gauge theories with N=3,4 and 6, normalized to 
their Stefan-Boltzman values, (b) The same for the entropy density. 



Other bulk thermodynamic quantities can be obtained from a knowledge of p{T) and A{T). 
Figure ^ shows a summary of our results for energy density, pressure and entropy density for SU(N) 
gauge theories, normalized to their Stefan-Boltzman limits. The Stefan-Boltzman values used are 
the results calculated for the free theory on the lattice, for the integral method [ [l3| ] : 
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where 



Sn^ 1 



(3.1) 
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As the figures show, even at temperatures close to 4 Tc, one reaches 85% or less of the Stefan- 
Boltzman value. The figures also show that, plotted as a function of temperature, the size of the 
1 /N^ corrections is small in all the three observables, except for very close to where 1 /N^ cor- 
rections ai^e visible in the energy density and entropy, as expected from the sizable 1 /N^ correction 
in the latent heat, Eq. (2.3). 



4. Summary and Discussions: 



We have presented results for bulk thermodynamic quantities for SU(N) gauge theories with 
N = 3, 4 and 6. Our results for latent heat of the transition and its dependence on N are summarized 
in Table |l| and Eq. (23). We find that the latent heat for the theory with three colors is substantially 
less than that with larger number of colors: the size of the leading correction in Eq. ( |2.3| ) is 
substantial for N=3. 

The bulk thermodynamic quantities are shown in Fig. |3[ It is found that even at temperatures 
close to 4 Tc, one reaches only ~ 85% of the Stefan-Boltzman value. We also see that except 
very close to Tc, the size of the leading 1/N^ correction is small, when we compai^e the bulk 
thermodynamic quantities measured at the same value of T /Tc. Interestingly, the scaling with the 
number of colors is better when looked at as function of T /Tc than when considered as function of 
the 't Hooft coupling. Figure ^ (a) shows the entropy density plotted against the 't Hooft coupling, 
defined through the V scheme and evaluated at the scale InT, where the running is done through 
our nonperturbative beta function. Considerable 1 /N^ correction is now seen, in particular for 
T < l.5Tc. Of course, the scaling with A will be different depending on the scheme chosen to 
define the coupling. The figure also shows s/sgg for the c/K = 4 supersymmetric theory [Q]. This 
theory is seen to have deviations from the Stefan-Boltzmann value of the same size as pure SU(N) 
theory at the highest temperatures, but a very different dependence on A. 
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Figure 4: (a) Entropy density, in units of the Stefan-Boltzman value, shown as a function of the 'tHooft 
coupling. The pink dotted line is the result for the JV — 4 supersymmetric SU(N) theory for large N [Q]. 
(b) p/psB shown as a function of e/Esb- The diagonal line is the line of conformal theories while, by 
construction, the Stefan-Boltzman limit is the point (1,1). The weak coupling results for the different theories 
are also shown. 

The large deviation of the thermodynamic quantities from the Stefan-Boltzmann value even 
at 4 Tc, and the fact that some strongly coupled conformal field theories show similar deviations 
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from Stefan-Boltzmann limit, have been used in the literature to speculate about a strongly cou- 
pled, conformal regime in pure SU(N) gauge theories. The thermodynamics results can be used to 
investigate the feasibility of such a phase. Following Ref. [14], in Fig. ^ (b) we plot the energy 



density vs. pressure, normalized by the corresponding Stefan-Boltzmann values. By construction 
the point at (1,1) is the Stefan-Boltzmann limit, while the diagonal line denotes conformality. Also 



shown are the weak coupling lines for the theories with the different number of colors QISQ. We 
find that for the temperature regime of interest to RHIC, ~ 2Tc, the theory is far from conformal. 
At higher temperatures, the theory approaches conformality, but it is closer to the weak coupling 
line than the conformal line, indicating the absence of a strongly coupled, conformal phase in the 
SU(N) gluon plasma. 

The computations were carried out on the workstation farm of the department of theoretical 
physics, TIFR. We thank Ajay Salve for technical support. 
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